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Abstract 

We study the entropy generation and particle production in scalar quantum field 
theory in expanding spacetimes with many-particle mixed initial states. The recently 
proposed coarse-grained entropy approach by Brandenberger et al. is applied to sys- 
tems which may have a non-zero initial entropy. We find that although the particle 
production is amplified as a result of boson statistics, the (coarse-grained) entropy 
generation is attenuated when initial particles are present. 
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One of the interesting features of quantized fields in a curved spacetime is that the con- 
cept of particles becomes very observer- dependent. For instance, in an expanding Universe 
spontaneous particle creation can occur. One defines generally a vacuum state such that all 
inertial observers in the past region agree that the spacetime looks empty of particles. As a 
result of the expansion of the Universe, the above vacuum state looks full of particles using 
modes natural to inertial observers in the far future region. Stated differently, a no-particle 
initial state can evolve to a many-particle state. However, since one starts with a pure state, 
one ends with a pure state. Thus there must be subtle correlations between the particles in 
the final state. In particular, there is no entropy production in this process even if lots of 
particles are produced. But, it may be that some of these subtle correlations are very diffi- 
cult to detect and/or that they may be quite sensitive to interactions between the produced 
particles. One may then contemplate a picture where such information about the system is 
considered to be "less relevant" and either discarded altogether or some kind of "statistical 
averaging" procedure is applied to it. This way one can try to associate a "coarse-grained" 
entropy to the final state of the system, hopefully in as natural way as possible. There has 
been a lot of work in this direction by Hu, Kandrup and collaborators 0. 

Recently, two different such approaches have been proposed. Brandenberger, Mukhanov 
and Prokopec (BMP) discussed in ||^, Q a coarse-graining procedure based on averaging 
over the so called squeeze angles which appear in the S-matrix of particle production. On 
the other hand, Gasperini and Giovannini together with Veneziano (GGV) P] related 
entropy generation to an increased dispersion of a superfluctuant operator. Both groups 
were especially interested in the entropy generation related to the production of gravitational 
waves and density fiuctuations in infiationary universe models. 

In this work, our aim is to follow the BMP approach and study the entropy generation 
starting not from an initial vacuum state with zero entropy, but allowing the system to be 
initially in some generic many-particle (mixed) state with non-zero entropy. If one starts with 
many bosons it is known 0] that the particle production will be amplified as a result of boson 
statistics, as one would expect. So, in general one can ask whether the entropy generation 
(in the coarse-grained sense) would also be amplified or not. Indeed, as a consistency check 
it is necessary to investigate if definitions of coarse-grained entropy will lead to a growing 
entropy even if initial state is allowed to be an arbitrary many-particle state with initial 
entropy. In ^ the GGV entropy was shown to be growing at least in certain classes of 
initial states. Interestingly, it was found that their entropy generation did not depend at all 
on the number of particles or entropy of the initial state. Here we will attempt to investigate 
the BMP entropy in similar situations. At least in the case of an initial density matrix where 
particles appear as pairs of opposite momenta, and initial entropy depends on the average 
occupation number per mode, we can show that the BMP entropy grows, though the entropy 
generation is attenuated. The BMP entropy does depend on the initial number of particles 
in a non-trivial way. In the end we comment briefiy on the case of an initial thermal density 
matrix. 

A scalar field in a D-dimensional curved spacetime is described by an action 

S = j d^xy^ - (m^ + ei?(x))02], (1) 

where R{x) is the Ricci scalar curvature of the metric and ^ is a coupling constant. Assume 
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that the metric depends exphcitly on time and that it is asymptotically flat in the far past 
and far future 

9fi,yix,t) C±rif,^ (2) 

as t — ^ ±00. In this case there are two natural ways to quantize the field in Heisenberg 
picture : 

cp = T.^i^f + h.c. (3) 

k 

(t^ = Y.^Tuf' + h.c. (4) 



where m^" looks like a plane wave in the far past, 

, t) oc e~H"*+^'^-"^ , t ^ -00 (5) 



and u°^^ looks like a plane wave in the far future, respectively 

[x,t) OC e k , t ^ +00 . (6) 

These solutions can be related via a Bogoliubov transformation, which can be given in terms 
of annihilation/creation operators as 

This transformation is generated by a S-matrix 



k k 



which has the explicit form 



(9) 

The factor , is the in-out vacuum amplitude. We use the convention of where 

the coefficients a have taken to be real. The S-matrix is known to generate a unitary 
transformation between the 'in' and 'out' representations if the gravitational field has a 
compact support For Robertson- Walker type universes the in-out vacuum amplitude is 



zero and the 'in' and 'out' representations are thus unitarily inequivalent. 
The S-matrix relates the in- and out-vacuum states in the following way 

I 0, in) = S I 0, out) 

^ exp{ha-'P]^,^,,af'''a^f^'} \ 0, out) . (10) 



detfa 
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This is the statement that an inertial observer in the far future region sees the in-vacuum 
state as full of out-particles. Similarly, the density matrix of the system expanded using 
in-modes (= p,) can be related to an expression using out-modes (= pf) as follows 

oo 

A = n H /(^fc) I n;:,in){in,nj: \ (11) 



k "fe 



nr=0 



Try . . . SaJ.°"*S-^S | O,out){out,0 \ S-^Saf'S"^ ■ ■ ■ Saf'S"^ 



k "fc ^fc- 



oo 



n fi^k)^ I n^, out) {out, I S ^ = Sp/S 



-1 



k 



nr=0 



Suppose now that the system is initially in an arbitrary many-particle state. In this state 
the average occupation number per mode (using in-modes) is given by 



n 



Trip^faf) . (12) 



In the far future region an inertial observer sees the average occupation number per mode 
using out-modes as 

— ^ T^^i'g^ a-l^'fout out 



n 



k TrpfTr[SpfS--al-^ar) (13) 

Using the cyclicity of the trace and the properties of the Bogoliubov transformation one can 
derive the relation between and to be 

k k 

This is the formula for "stimulated emission" It tells us that even if the spontaneous 
creation of particles is weak, | /J^^. 1, the particle production can become arbitrarily 
large, 

An,^ ^ n{ - nj: =| P^,^. \' (1 + 2n» ^ large (15) 

if the initial average occupation number per mode is arbitrarily large. This amplification 
of particle production is a result of the boson statistics of the particles. For fermions the 
particle production would be attenuated |]lT] . 

Let us now discuss for notational simplicity metrics of the form 

ds'^ = dt^ - a^{t)£x^ , (16) 

where a{t) is a scale parameter of the universe. We again just require that a{t) — ^ a± 
asymptotically as t ^ ±oo. As a result of the invariance under spatial translations, the 
Bogoliubov coefficents can be written as 

ar- = «r5r cosh rr^r - 

kp k k,p k k,p /I y^j 

Pkp = Pkk-p^^'''^^^^r^kh,-p- ^ ' 

The parameters rg, Oj: are called squeeze parameter and squeeze angle in the Quantum Optics 
litterature, and the S-matrix is called a two- mode squeeze operator ||12|. If one starts with 
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a vacuum state, the final state (10) is called a squeezed vacuum. In the initial vacuum case, 
if one expands the corresponding Sp/S~^ in the 'out' basis of energy eigenstates, one finds 
that the off-diagonal components of Sp/S~^ have an oscillatory dependence of the angles 
Ok- In the BMP coarse-grained entropy approach it is assumed that these angles represent 
irrelevant information about the system {e.g., in the sense that they would be very difficult 
to measure) and they are therefore averaged over. After the averaging only the diagonal 
elements of Sp/S^^ survive and one then defines a coarse-grained entropy with the resulting 
reduced density matrix p^ed with the usual formula S = —kBTr{pred^ogpred)- The result is 



^ffi ^ ^ gf^fi ^ f^^ ^(cosh^ log cosh^ — sinh^ log sinh^ r£) 
k ^ 

- ks E[(^f + 1) log(^f + 1) - (^f ) log(^f )] ' 



where the notation nl'^ means the LHS of (14) in the case of an initial vacuum state. Now we 
turn to consider initial density matrices pi which can describe generic many-particle states 
with non-zero entropy. Let us assume that the initial density matrix has the form 

oo 

pi= n E/fcH I ^fc'"-fc'^")(^'^'"'fc''^-fc I ' (19) 

k,{k^>0) 

where the coefficients /(^g) are of the form f{n^ = (n*-)"/(n*- + 1)"^"^. That is, we start 
with a many-particle state where particles appear in pairs of opposite momenta, with an 
initial average occupation number spectrum = and with (ordinary) entropy given by 
— kBTr{pi log Pi) . Writing the initial entropy in more explicit form, it is 

S'^ H sl = kB E [(rii: + l)log(n^g + l)-(n|)log(n|)] . (20) 

fc,(fcz>0) k,{k^>0) 

When we expand the resulting final density matrix Sp/S~^ (11) in 'out' energy eigenstates, 
we find that also in this case the off-diagonal elements have an oscillatory dependence on 
the angles Oj:. Therefore, following the BMP approach and averaging over the angles, only 
the diagonal elements will survive. Thus the reduced density matrix of the final state has 
the form 

oo 

Pred= n E /fcH I "-fc' ^-fc' O"^) ^fc' ^-fc I ' (21) 

k,{k^>0) "=0 

where fj^{n) = {nj:,n_j:,out \ SpfS~^ \ out,nj:,n_j:) . After some effort, one can show that 
the coefficients have the form 

= rL un+i ' (22) 
where is the LHS of (14) with the of (19). Thus, the final coarse-grained entropy is 

E 4 = kB E [(r^^ + 1) log(ni + 1) - (n^ log(nf.)] . (23) 

k,{k^>0) k,(k^>0) 
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The entropy depends only on the occupation number spectrum of particles in the final state. 
This result is in agreement with a similar formula given in by a more heuristic argument 
to define entropy of a statistical system with a definite spectrum which is valid both in and 
far out of thermodynamical equilibrium. Let us now compare the entropy generation per 
mode in the initial vacuum and initial many-particle cases. Denote 

^os^k = 4'' - (24) 
= 4 - 4 ' (25) 

where (24) applies to the former case and (25) to the latter case. As a first consistency check, 
we find that 

As^, > (26) 

so that the coarse-graining led to a growing entropy in our many-particle case. However, as 
we compare (25) and (24) we find that 

As^ < Aos-^ ; (27) 

i.e., the entropy generation is attenuated if one starts with many particles present in the 
mode k. The equality holds iff = 0. The difference As^ — AgSg depends symmetrically 
on n*. and | /?g P and decreases monotonically as either variable increases. Further, as both 
variables approach infinity, 

Asg- AqS^ ^ log2- 1 -0.31 (28) 

asymptotically. This finite value is the maximum difference between the generated entropies 
per mode. Thus, unlike the GGV entropy, the BMP entropy is not independent of the 
number of particles in the initial state, but has some 'memory' about the initial occupation 
numbers. Since entropy is a measure of loss of information, it would appear that more 
information about the initial state of the system is conveyed to the final coarse-grained state 
when stimulated emission dominates the spontaneous particle production (since entropy 
generation is attenuated). 

Next case to be investigated would be an initial thermal density matrix 



Pi = n^fc ^^^p(-/^' 

it 

This situation is somewhat trickier to deal with, for the following reason. Initially, the 
particles of opposite momenta are uncorrelated. However, in the expansion of the universe 
the particles are produced in pairs of opposite momenta. This induces correlations between 
the opposite momenta in the final density matrix. It would have the form 

oo 

SpfS'^ = n E 4(^' I ^fc' ^-h ^"^) (^"^' ^'fc' ^'-k I ' (29) 

k,{k^>0) n,m,n',m'=0 

where 

f^{n,m;n' ,m') = 



ZkZ., 
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Again, one can see that the n ^ n' oi m ^ m' components have an oscillatory dependence of 
the squeeze angles and they vanish in the coarse graining. However, the diagonal coefficients 
(those of the reduced density matrix) will have a form /^(n, m) where the dependence on n 
and m does not factorize. Hence the opposite momenta have aquired correlations through 
the particle production and the reduced density matrix is not of the same type as the initial 
one. As advocated in [^], one would like to ignore the correlations between different modes. 
In [0 this problem was avoided by replacing the two-mode squeeze operator by a one-mode 
squeeze operator, but this does not appear entirely satisfactory, since the particles are then 
not created in the correct way as pairs of opposite momenta. We would like to propose 
that the correlations between opposite momenta could be ignored by proceeding to define 
f^iji) = J2m fki^j^) = I]n /fcl'^) '"^)- Then we would define the final reduced 

density matrix to be 

oo 

Pred= n Mn)f_^{m)\n^,m_^,out){out,nj:,m_^\ , (31) 

fc,(A:^>0) n,m=0 

which is of the same type as the initial density matrix. Now the final entropy would be given 
by 

= 'kBY.T.Un)\ogUn) . (32) 

k " 

Unfortunately, at the present we do not have explicit formulas for the coefficients fj:{n) or 
the final entropy. It would be very interesting to see if the resulting expressions could depend 
on the final average occupation number in the same fashion as in the earlier case. We hope 
to be able to return to this question in the future. 

Finally, let us clarify that even if we found a different result as in the GGV approach, 
namely an entropy that depends on the number of particles in the initial state, we are not 
arguing that it would mean that the BMP approach is 'better' than the GGV approach. 
As stated in [0], it is good to have different definitions of entropy, corresponding to loss of 
different information about the system. Both BMP and GGV approaches have the virtue 
of giving the correct average occupation number for particles in the final state. Otherwise 
the GGV approach appears to discard information about the system a bit more generously, 
since it leads to a greater growth of entropy. 
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